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1.  Introduction 


The  familiar  intensity  distribution  associated  with  Fraunhofer  diffraction  from  a  slit  or  aperture 
may  be  observed  not  only  in  the  far  field,  but  also  in  the  focal  plane  of  a  well-corrected  lens. 
Accordingly,  an  accurate  determination  of  damage  thresholds  for  focal  plane  elements  must 
correctly  account  for  diffractive  effects. 

Of  course,  a  reported  value  of  the  fluence  (incident  energy  per  unit  area)  at  which  an  optical 
component  experiences  damage  is  truly  meaningful  only  if  that  fluence  is  appropriately  defined. 
Depending  on,  for  instance,  the  size  or  thermal  properties  of  the  component  in  question,  the  peak 
fluence  to  which  the  component  is  subject  may  prove  a  less  useful  measure  of  the  damage 
threshold  than  the  average  fluence  impinging  on  the  component  or  on  some  portion  thereof. 
Thus,  in  cases  in  which  an  average  value  of  the  fluence  is  deemed  to  provide  the  most 
meaningful  measure  of  the  damage  threshold,  one  must  exercise  a  certain  amount  of  care  in  the 
choice  of  convention  regarding  the  area  over  which  average  fluences  are  to  be  calculated. 

In  this  technical  note,  we  review  the  physics  of  Fraunhofer  diffraction  from  a  uniformly 
illuminated  circular  aperture  (a  top  hat  beam)  and  describe  the  computation  of  peak  and  average 
fluences  in  the  presence  of  diffractive  effects. 


2.  Diffraction  From  a  Circular  Aperture 


We  consider  the  situation  illustrated  in  Figure  1:  a  lens  of  focal  length  L  situated  in  a  circular 
aperture  of  radius  a ,  upon  which  a  planar  wave  front  is  incident.  (In  practice,  the  incident 
radiation  is  provided  by  a  high-quality  Gaussian  laser  beam  that  overfills  the  aperture  to  the 
extent  that  the  wave  front  can  be  well  approximated  as  planar.)  The  incident  radiation  is 
monochromatic  with  wavelength  A;  k  =  2n / A  is  the  wavenumber.  Let  p  denote  the  radial 

coordinate  in  the  image  plane  and  for  later  convenience  define  w  =  sin#  =  p! -Jp2  +  L 2  . 


The  irradiance  profile  in  the  image  plane  is  given  by  the  familiar  result  of  Airy,1 


/(sin  6)  =  1(0) 


2J)  (for  sin  0) 
ka  sin  0 


indicated  by  the  black  line  in  Figure  2.  Dark  rings  are  observed  when  ka  sin  <9,  the  argument  of 
the  Bessel  function  in  equation  1,  assumes  a  value  associated  with  a  root,  i.e.,  when  ka  sin#= 
3.832,  7.016,  10.174, ...,  or  equivalently,  when  sin#=  0.610  A/a,  1.116  A/a,  1.619  A/a,  etc.  For 
the  innermost  dark  ring,  0  « 1 ,  so  sin  6  *  tan  0  =  p! L ,  and  one  can  write  the  following 


'Airy,  G.  B.  Tram.  Camb.  Phil.  Soc.  1835.  5.  283. 
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Figure  1.  Diffraction  from  a  circular  aperture. 


Figure  2.  The  Airy  pattern  (with  “best  fit”  Gaussian). 

well-known  expression  for  the  radius  of  the  first  dark  ring  in  terms  of  the  ratio  L/2a ,  the  /-number 
(//#)  of  the  system:  px  «  0.610/1—  =  1.22  X{f  /#). 

Using  a  very  general  argument  involving  energy  conservation,  one  can  show  that  the  on-axis 
irradiance  7(0)  is,  to  a  very  good  approximation,  given  by  the  following  simple  formula: 


2 


(2) 


7(0)  = 


m2  P 

¥l2  ’ 


in  which  P  is  the  total  power  passing  through  the  aperture.  Equation  2  is  obtained  from  an  ex¬ 
pression  for  the  irradiance  that  is  derived  with  the  restriction  that  sin  6  « 1 .  This  expression  is 
integrated  over  the  image  plane.  However,  the  error  made  in  extending  the  integration  over  the 
entire  image  plane  is  negligible,  since  the  integral  is  dominated  by  the  contribution  from  the 
paraxial  region,  in  which  the  restriction  holds;  at  large  angles  6  from  the  optic  axis,  the 
irradiance  is  vanishingly  small.2 

Integration  of  equation  1  over  all  time  yields  the  total  fluence  F(sin#) ,  which  is  simply  the 
product  of  the  irradiance  (equation  1),  the  pulse  width,  and  a  dimensionless  constant  dependent 
on  the  pulse  shape.  For  instance,  for  a  single  pulse  that  is  Gaussian  in  time  with  HW(l/e2)M 

pulse  width  r ,  the  fluence  is  given  by  /'’(sin  6)  =  4ft  z  /(sin  6) ,  while  for  a  square  pulse  of 
duration  T ,  it  is  simply  T  /(sin  6) .  The  spatial  profile  of  the  fluence  is  identical  to  that  of  the 
irradance: 


F(sin  0)  =  F(  0) 


2Jx(ka  sin  6) 
kasmd 


(3) 


and  F( 0) ,  the  fluence  at  the  focal  point,  is 


F(  0)  = 


7ta2E 
A2  L2 


7ZE 

4/1 2(//#)2  ' 


(4) 


In  equation  4,  E  represents  the  total  energy  passing  through  the  aperture  and  /  /#  is  the  /-number 
of  the  system,  l/la . 

The  appearance  of  A2  in  the  denominator  of  equations  2  and  4  might  at  first  seem  puzzling.  To 
make  the  presence  of  this  factor  more  plausible,  one  might  recall  that  the  effects  of  diffraction 
from  the  aperture  diminish  in  importance  at  short  wavelengths.  Thus,  for  small  A ,  one  would 
expect  to  see  an  increase  in  the  amount  of  energy  that  the  lens  concentrates  at  focus,  just  as 
equation  4  predicts. 


3.  Encircled  Energy  and  Average  Fluence 


In  the  determination  of  damage  thresholds  for  certain  types  of  focal  plane  elements,  the  peak 
fluence  F(0)  may  prove  less  useful  than  the  average  fluence  incident  on  a  given  region.  This  is 
defined  simply  as  the  ratio  of  the  total  energy  deposited  in  the  region  to  the  area  of  the  region.  In 
the  present  context,  in  which  the  setup  is  cylindrically  symmetric  about  the  optic  axis,  it  is 
natural  to  consider  a  disk-shaped  region  of  the  image  plane  with  the  focal  point  at  the  center. 


2 

Bom,  M.;  Wolf.  E.  Principles  of  Optics,  6th  ed.;  Cambridge  University  Press:  Cambridge,  UK;  1980;  p  386. 
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The  average  fluence  incident  on  a  disk  of  radius  p0  (which,  seen  from  a  point  at  the  center  of  the 
aperture,  subtends  the  angular  radius  w0  =  pQ  /  -\] p0  +  L  )  is 


F(w0)  = 


E{w  q) 

2  ? 


(5) 


where  E{w0 ) ,  the  total  energy  deposited  in  the  disk,  is  the  encircled  energy  in  the  solid  angle 
subtended  by  the  disk,  as  seen  from  the  center  of  the  aperture.  E{w0 )  is  the  integral  of  the 
fluence  over  the  disk. 


A) 

E{w0 )  =  In  ^F{w)pdp  =  2^F(0) 
0 


I2  "(r  dw  [2 J,  {kaw)]2 

k2a2  jl-w2  w 


(6) 


Since  w2  <  1 ,  the  factor  [1  -  w2  ]_1  in  the  integrand  of  the  latter  equality  can  be  replaced  by  the 
geometric  series  Iw2" .  Doing  so,  and  substituting  for  F{ 0)  from  equation  4,  one  obtains 


oo  w0 

E(w0)  =  ±E^]cbvw2n-l[2J1(kcnv)]2 

n= 0  0 

=  2E  {  \dw  —  ^caw\  +  \dw  wJ2 ( kaw )  +  . . . 

[  0  w  0 


(7) 


The  integrations  may  be  performed  by  making  use  of  the  recursion  relations  for  the  Bessel 
functions.  The  result,  expressed  as  the  fraction  )  of  the  total  incident  energy  contained 
within  a  circle  in  the  image  plane  that  is  centered  at  the  focal  point  and,  seen  from  the  center  of 
the  aperture,  subtends  the  solid  angle  2/rw0 ,  is 


=  1  -Jo  (kanv 0 )  -  J2  {kaw  0 )  (g) 

+  wl  \j\  (*oiv0  )-J0  ( kaw0 )  J 2  {kaw0 )] + o(w„  [kaw^1 ). 

To  the  best  of  our  knowledge,  this  is  a  new  result,  but  it  is  of  little  practical  interest,  since  the 
version  of  this  formula  published  by  Lord  Rayleigh  in  1881,3 

s(w0  )  a  1  -Jo  ( kaw0 )  -  J2  {kaw 0  ) ,  (9) 

though  valid  only  to  zeroth  order  in  w0 ,  suffices  completely  for  most  applications.  It  is  plotted 
in  Figure  3. 


3  Rayleigh.  J.  W.  S.  Phil.  Mag.  1881, 11. 214. 
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Figure  3.  Fraction  of  total  energy  contained  in  a  disk  in  the  focal  plane  that  is 
centered  at  focus  and  subtends  solid  angle  2n  sin  #at  the  center  of  the 
aperture. 


Written  in  terms  of  £(w0) ,  the  average  fluence  is 


F(w0)  = 


Es^w q) 
XPo 


=  Es(w0) 


i-W 

7Vw\L 2 


Es{w  q) 
nw\l}  ’ 


(10) 


where  the  last  expression  is  valid  for  0  «  1 .  Here,  the  disk  radius  pQ  (or  equivalently,  w0 ,  the 

angular  radius  subtended  by  the  disk)  is  arbitrary,  and  various  workers  may  make  different 
choices  for  the  size  of  their  “standard  disk.”  A  particularly  natural  choice  is  to  take  the  disk  to  be 
the  central  bright  disk  of  the  Airy  pattern,  the  region  bounded  by  the  first  dark  ring  at 
w0  =0.610/1  la.  The  central  Airy  disk  contains  84%  of  the  energy  passing  through  the  aperture. 

Alternately,  one  might  choose  iv0  to  lie  at  the  half-width  at  half-maximum  (HWHM)  of  the 

central  Airy  peak  or  the  half-width  at  1  /  e2  of  maximum  (HW(  1  /  e2  )M)  of  the  central  peak,  this 
later  convention  being  justified  by  the  fact  that  although  the  focal  plane  fluence  profile  (equation 
3)  is  not  Gaussian,  the  central  peak  is  reasonably  well  approximated  by  a  Gaussian,  as  Figure  2 
illustrates.  A  disk  of  radius  equal  to  the  HWHM  of  the  central  peak  has  angular  radius 
w0  =  0.257  X  la  and  contains  47%  of  the  energy  passing  through  the  aperture,  while  a  disk  of 

radius  equal  to  the  HW(  1  /  e2  )M  of  the  central  peak  has  iv0  =  0.4 1 1 X  la  and  contains  77%  of  the 
energy. 


4.  Numerical  Example 


To  better  illustrate  the  extent  to  which  the  numerical  value  of  a  reported  fluence  depends  on 
particular  convention  used,  we  consider  an  aperture  10  cm  in  diameter  admitting  1  p  J  of 
radiation  at  532  nm.  For  an  JI5-,  an  f/2A-,  and  an  fl\. 2-system,  we  compute  the  following: 

•  F(0) ,  the  fluence  at  the  geometrical  focal  point; 
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•  F(w0  =  0.6 1 0  X  /  a) ,  the  average  fluence  over  the  central  disk  of  the  Airy  pattern;  and 

•  F(w0  =  0.41  l/l /a) ,  the  average  fluence  over  a  disk  of  radius  equal  to  the  HW(  1  /e2)M  of 
the  central  Airy  peak. 

The  results  are  summarized  in  Table  1. 


Table  1.  Fluence  values  for  a  10-cm  aperture  admitting 
1  |iJ  of  energy  at  A  =  532  nm. 


Fluences 

(J/cmz) 

//  5 

ft  2.4 

// 1.2 

Fluence  at  geometric  focal  point 

11.1 

48.2 

193 

Average  fluence  of  central  Airy 
disk 

2.53 

11.0 

44.0 

Average  fluence  over  a  disk  of 
radius  the  HW(  l/e2)M  of  the 
central  airy  peak 

5.10 

22.1 

88.5 

5.  Practical  Formulae 


We  conclude  with  a  compilation  of  practical  formulae  assembled  for  the  benefit  of  workers 
measuring  damage  threshold  fluences  for  focal  plane  elements.  The  formulae  for  average  flu¬ 
ence  are  obtained  by  applying  the  general  results  equations  9  and  10  to  particular  choices  for  a 
“standard  disk”  discussed  at  the  end  of  section  3.  As  before,  E  represents  the  total  energy 
passing  through  the  aperture,  (//#)  =  LAa  is  the/-number  of  the  system,  and  p0  is  the  radius  of 

the  standard  disk  in  question. 

•  From  equation  4,  the  peak  fluence,  i.e.,  the  fluence  at  the  geometric  focal  point: 

m=0JS5J(7W 

•  Average  fluence  over  the  central  disk  of  the  Airy  pattern: 

E  E 

F  (central  disk)  =  0.267— j  =  0.179  vf- 

Po  4  \J  '#) 

•  Average  fluence  over  a  disk  of  radius  equal  to  the  HWHM  of  the  Airy  peak: 

F(HWHMdisk)  =  0.15lA  =  0.572-^. 
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•  Average  fluence  over  a  disk  of  radius  equal  to  the  HW(  1  /  e2  )M  of  the  Airy  peak: 

F  (HW(l/e2  )M  disk)  =  0.244^-  =  0.361 — — — — — . 

Po  A2  {f/#Y 
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